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Abstract 

Afforestation is a mitigation option to reduce the increased atmospheric carbon dioxide levels as well as the predicted high 
possibility of climate change. In this paper, vegetation survey data, statistical database, National Forest Resource Inventory 
database, and allometric equations were used to estimate carbon density (carbon mass per hectare) and carbon storage, 
and identify the size and spatial distribution of forest carbon sinks in plantation ecosystems in sand source areas of north 
Beijing, China. From 2001 to the end of 201 0, the forest areas increased more than 2.3 million ha, and total carbon storage in 
forest ecosystems was 1 73.02 Tg C, of which 82.80 percent was contained in soil in the top 0-1 00 cm layer. Younger forests 
have a large potential for enhancing carbon sequestration in terrestrial ecosystems than older ones. Regarding future 
afforestation efforts, it will be more effective to increase forest area and vegetation carbon density through selection of 
appropriate tree species and stand structure according to local climate and soil conditions, and application of proper forest 
management including land-shaping, artificial tending and fencing plantations. It would be also important to protect the 
organic carbon in surface soils during forest management. 
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Introduction 

Terrestrial vegetation as an important component of global 
carbon cycle stores more than 600 Gt C and annually exchanges 
about 1 0 percent of the carbon storage with the atmosphere via 
photosynthesis and respiration [1]. Among all vegetation types, 
forests play an important role in regional and global carbon cycles 
because they store large quantities of carbon in vegetation and 
soils; they become atmospheric carbon sources following human or 
natural disturbance, and become carbon sinks during re-growth 
after disturbance [2]. Forestry projects improve tree coverage on 
lands that have not supported forest growth in historical time 
(afforestation), through direct seeding and planting, and through 
plantation forestry [3]. Enhancing carbon sequestration by 
increasing the plantation area is regarded as a potential measure 
to mitigate atmospheric CO2 concentration rise and prevent 
global warming [4] . The effects of afforestation on total ecosystem 
carbon storages have been widely studied in recent years [5-1 1], 
and consequently, global, national, and regional carbon densities 
databases in various above- and belowground forest pools [12-14] 
in different forest types [15], land uses [16-17], and forest areas 
[18-20] have been established. National Forest Resource Inven- 
tory (NFRI) database (2004-2008) from the Forestry Ministry of P. 
R. China shows that 80 percent of plantations in China are young, 
middle, and near-mature forests [14,20]. To understand the effects 
of afforestation projects especially at early implementation stage on 



carbon dynamics, estimation of carbon storages in young 
plantation ecosystems is urgently needed. 

In many areas of North China, the intensified spring wind 
possibly due to climate changes, and the inappropriate human 
activities including relentless land cultivation, deforestation, and 
overgrazing, have led to a continued loss of vegetation cover and 
topsoil [21-22]. Facing such serious ecological problems, the 
Chinese government has invested huge amounts of money into 
Key Forestry Programs including Grain for Green Project, Three- 
North Shelter Forest System Project, and Beijing-Tianjin Sand 
Source Control Project (BTSSCP) [23]. The implementation of 
these programs has established 62 million ha of afforestation by 
the end of 2010, with a total forest cover of 20.36 percent. 
Plantations can improve vegetation cover and carbon sequestra- 
tion, and weaken the transport of wind-born dust and sand to 
other countries by controlling soil erosion [24—25]. In 2000, a 
national ecological construction program, the BTSSCP, was 
launched to improve the ecological environment in Beijing, 
Tianjin, and the surrounding areas and thereby to prevent the 
hazards of sand storms by forest recovery and protection [26]. 
Estimating carbon storages in plantation ecosystems in sand source 
areas plays an important role in objectively evaluating the eco- 
environmental effects of treatment projects and affirms the 
contribution of Chinese terrestrial ecosystems to the global carbon 
cycle. 

In this study, vegetation survey database, allometric equations, 
statistical database of annual planting species in sand source areas 
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from 2001-2010, and 2004-2008 NFRI were used to determine 
(1) carbon storages and carbon densities and (2) the size and spatial 
distribution of carbon sinks over the past decade in plantation 
ecosystems in sand source areas of north Beijing. 

Materials and Methods 

Site description 

The study area (38°51'-45°11' N, 109°53'-120°58' E) is 
located in the sand source areas (29.6 million ha) of north Beijing 
(Fig. 1). It covers Haihe River Plain, Taihang Mountain, Yanshan 
Mountain, and the Inner Mongolian Plateau. According to the 
distribution laws of bioclimatic zones and geomorphic types at a 
regional scale, the study area was divided into water resource 
protection zones in Yanshan mountainous and hilly region 
(WRPZYMHR), desertificated land zone in agro-pasture region 
(DLZAPR), and Otingdag sandy land zone (OSLZ) (Table 1). The 
elevation ranges from 10 to 2000 m a.s.l. The climate includes 
warm temperate semi-humid, temperate semi-humid, temperate 
semi-arid, and temperate arid zones. Annual rainfall averages 
459.5 mm, with 500-600 mm in WRPZYMHR, 250-450 mm in 
DLZAPR, and 300^00 mm in OSLZ. Annual potential evapo- 
transpiration averages 2110 mm with 1600-2200 mm in 
WRPZYMHR, 1800-2500 mm in DLZAPR, and 2000- 
2500 mm in OSLZ. Annual temperature averages 7.5°C with 
4°C in WRPZYMHR, 2-8°C in DLZAPR, and 1.6-7.0°C in 
OSLZ. The soils were classified by Chinese Soil Classification as 
chernozems, castanozems, brown earths, calcareous soil, and 
lithosols [27], corresponding to mollisols, calciustoll, alfisols, 
andisols, and entisols in United States Department of Agriculture 
(USDA) Soil Taxonomy [28]. 

There are more than 2300 plant species belonging to 130 
families and 670 genera. The major families include Compositae, 
Gramineae, Ranunculaceae, Cyperaceae, and Rosaceae. The 
natural vegetation on Taihang and Yanshan Mountains is 
dominated by temperate and warm temperate deciduous broad- 
leaved forests such as deciduous oak forests and small leaf 
deciduous forests. The dominant tree species include Quercus 
wutaishanica, Q. mongolica, Q. aliena, Q. acutissima, Q. variabilis, Betula 
platyphylla, Populus davidiana, and Ulmus pumila. However, the 
majority of existing vegetations are P. davidiana and B. platyphylla 



forests, Vitex negundo var. heterophylla, Lespedeza bicolor, and Armeniaca 
sibirica deciduous shrubs. The plantations are dominated by Pinus 
tabulaeformis and Larix gmelinii. The natural vegetation on the Inner 
Mongolian Plateau is dominated by shrub-herb communities, 
including Stipa grandis, S. hylovii, S. breviflora, and Artemisia frigida. 
The plantations only account for a smaller proportion of forests 
and are dominated by broad-leaved forests and xerophilous 
shrubs. The planted tree species in the sand source areas are listed 
in Table 1. 

Field measurements and chemical analysis 

After reconnaissance surveys with permission from County 
Forestry Bureaus of Hebei province, Shanxi province, and Inner 
Mongolia Autonomous Region, the communities with broad- 
leaved forests (dominated by P. davidiana, B. platyphylla, and U. 
pumila), coniferous forests (dominated by P. tabulaeformis, L. gmelinii, 
Pinus sylvestris var. mongolica, and Picea asperata), economic forests 
(dominated by A. sibirica and Malus pumila), and shrubs (dominated 
by Hippophae rhamnoides, Caragana korshinskii, Salix gordejevii, and 
Hedysarum laeve) were selected as the representative forest types in 
WRPZYMHR, DLZAPR, and OSLZ. The barren hiUs surround- 
ing the main vegetation types were selected to represent the state 
before afforestation. Four sample plots of 20x20 m each were 
systematically laid out on each stand (13 forest types x4 replicates 
x3 regions = 156 in total). Four sample plots of different ages (3-, 
5-, 8-, and 10-year-old) on each forest type (13 forest types x4 
stand age = 52 in total) were selected. Stand age was not 
replicated in our vegetation survey because it was difficult to find 
replicate stands of the same age with similar species composition, 
soil conditions, and environmental conditions in two regions [29]. 
Physiographic factors including longitude, latitude, elevation, 
slope, and slope aspect across forest types were measured by 
GPS (UniStrong G330) and a clinometer. The species of all trees 
falling within a sample plot were identified according to Flora 
Reipublicae Popularis Sinicae [30]. The abundance, height and 
cover of each species, the number of individuals, the breast height 
diameter of broadleaved and coniferous trees, and the basal 
diameter of economic trees and shrubs were measured. Biomass of 
trees was determined using destructive sampling techniques (i.e. 
total harvest including roots). Six standard trees of each species at 
each stand age were harvested for measurement of above- and 




Figure 1. The location of plantation ecosystems in sand source areas of north Beijing, China. I: Water resources protection zone in 
Yanshan mountainous and hilly region; II: Desertificated land zone in agro-pasture region; III: Otingdag sandy land zone. 
doi:1 0.1 371 /journal.pone.0082208.g001 
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Table 1. The planted tree species in sand source areas of North Beijing, China. 





Vegetation type 
region 


Longitude and latitude 


Area 

(million ha) 


Planted tree species 


Water resources 
protection zone in 
Yanshan 

mountainous and 
hilly region 


39°34'53"-42°37 
119014' 5"E 


43"N 113°54'21"- 


7.5 


Populus hopeiensis, Robinia pseudoacacia, Ulmus glaucescens, Castanea mollissima, 
Betula platyphylla, Salix matsudana var. matsudana, Quercus mongolica, Salix 
babylonica, Acer mono, Platycladus orientalis, Pinus sylvestris var. mongolica, Pinus 
tabulaeformis, Larix gmelinii, Picea asperata, Armeniaca sibirica, Ziziphus jujuba, 
Malus pumila, Pyrus hopeiensis, Crataegus pinnatifida, Lycium chinense, Caragana 
korshinskii, Hippophae rhamnoides, Astragalus adsurgens, Amorpha fruticosa, Tamarix 
chinensis, Lonicera Japonica, Rhus typhina, Vitis vinifera, Vitex negundo var. 
heterophylla, Elaeagnus angustifolia, Salix cheilophila 


Desertificated land 
zone in agro-pasture 
region 


38^51 '26"-42 D 21 
1164' 10"E 


29"N 109*53 '4"- 


6.4 


Populus davidiana, Salix matsudana var. matsudana, Robinia pseudoacacia, Ailanthus 
altissima, Ulmus glaucescens, Juglans regia, Syzygium aromaticum, Larix gmelinii, 
Pinus sylvestris var. mongolica, Platycladus orientalis, Pinus tabulaeformis, Amygdalus 
davidiana, Armeniaca sibirica, Ziziphus jujuba, Pyrus hopeiensis, Prunus salicina, Vitis 
vinifera, Amorpha fruticosa, Hippophae rhamnoides, Spiraea salicifolia, Lycium 
chinense, Tamarix chinensis, Salix cheilophila, Lonicera Japonica, Medicago sativa, 
Ostryopsis davidiana, Caragana korshinskii 


Otingdag sandy land 
zone 


41°17'17"-45 c 11 
120°58'40"E 


23"N 112°36'14"- 


15.7 


Populus davidiana, Populus alba var. pyramidalis, Robinia pseudoacacia, Salix 
matsudana var. matsudana, Salix matsudana, Amorpha fruticosa, Acer mono, Acer 
truncatum, Ulmus glaucescens, Morus alba, Betula platyphylla, Xanthoceras sorbifolia, 
Pinus tabulaeformis, Platycladus orientalis, Picea asperata, Larix gmelinii, Pinus 
sylvestris var. mongolica, Armeniaca sibirica, Malus pumila, Ziziphus jujuba, Malus 
asiatica, Pyrus hopeiensis, Amygdalus davidiana, Prunus salicina, Ziziphus jujuba var. 
spinosa, Lycium chinense, Elaeagnus angustifolia, Caragana korshinskii, Hedysarum 
fruticosum var. mongolicum, Caragana microphylla, Hippophae rhamnoides, Vitis 
vinifera, Salix gordejevii, Astragalus adsurgens, Tamarix chinensis, Hedysarum laeve 



doi:1 0.1 371 /journal.pone.0082208.t001 



belowground biomass and tissue carbon content. The above- and 
belowground biomass was separated by tissue types (e.g. stem, 
branch, foliage, and root). The trees were cut at 20 cm above- 
ground. Prior to branch removal, four to six branches from bottom 
to top along the crown were sampled from each tree, and foliage was 
collected from each branch. All branches were then clipped off the 
tree. The stem of each tree was cut into 2-m-long sections using a 
chainsaw. The entire roots of the sampled trees were excavated and 
sub-samples were used to determine total root biomass. Fresh 
organs were weighed in the field; 500-1000 g of fresh samples were 
placed into nylon bags, and refrigerated and transported to our 
laboratory. Understory biomass was also measured using destructive 
sampling techniques (i.e. total harvest including roots), within three 
2x2 m, three lxl m, and three 20x20 cm randomly selected 
microplots which were established within each 20x20 m sample 
plot for data collection from shrubs, herbs, and litters. In each shrub 
or herb microplot, whole ground vegetation including roots (branch, 
foliage, and root for shrubs; above- and belowground parts for 
herbs) was harvested. Forest floor components (coarse wood, litter, 
and fragmentation layer) were collected within these 20x20 cm 
microplots. All samples of shrubs, herbs, and litters were weighed in 
the field, placed into nylon bags, and transported to the laboratory. 
All samples were dried in a forced-air oven at 65°C for 24 h. After 
measurement of dry weight, each sample was ground to pass a 1 00- 
mesh screen. Plant organic carbon was measured by the 
K 2 Cr 2 0 7 +H 2 S0 4 digestion method [31]. 

Soils were sampled at five depths (0-10, 10-20, 20-40, 40-60, 
and 60-100 cm) on five random locations along a diagonal 
transect within each of the 20x20 m plots for measurement of 
carbon content (2.0 kg) and soil bulk density (100 cm' per sample). 
For measurement of soil bulk density along with soil profile, the 
soil samples were weighed immediately and transported to the 
laboratory, where they were oven-dried at 105°C for 24 h and 
reweighed for calculation of soil bulk density and soil water 
content. To measure carbon content, the soil samples were pooled 



by the sampling location and layer, air-dried and ground to pass a 
100-mesh screen prior to analysis. Soil organic carbon (SOC) was 
also measured by the K 2 Cr 2 07+H 2 S04 digestion method [31]. 

Biomass estimation 

In this paper, forest carbon storage CTj (Mg ha ') for 
broadleaved, coniferous, and economic forests and shrubs was 
calculated separately as [32]: 

2 2 
CT f = £ (C n B„) + Y, {CujBuj) + 
t=i j=\ 

Y, {C Fm B Fm ) + (W 2 C s „BD Sn d) 

in — 1 n — 1 

where i is the tree tissue type (i.e. the above- and belowground 
parts), Ct (%) and B T (Mg ha~ ) are carbon content and biomass 
of tree tissues respectively; j is the understory component (i.e. the 
above- and belowground parts), Cy (%) and B L r (Mg ha ) are 
carbon content and biomass of understory components respec- 
tively; m is the forest floor component (i.e. coarse wood, litter, and 
the fragmentation layer), Cp (%) and Bp (Mg ha ) are the carbon 
content and the biomass of the forest floor component respectively; 
n is the layer of mineral soil, C$ (%) and BD S (g cm -3 ) are the 
carbon content and the bulk density of the measured soil layer 
respectively; d is the depth of the measured soil layer. 

The relationships between timber volume, tree carbon density and 
stand age were used to estimate above- and belowground biomass 
and carbon stocks for trees. Radial growths along the longest, 
shortest, and intermediate radii on each section were determined to 
calculate the stem volume over the bark [33]. NFRI database 
provides an average stand volume for different site conditions or 
forest management practices, and provides the stand volumes for 
some tree species or species groups of plantations [34]. Based on 
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Coniferous forests 




4 6 
Time (year) 

Figure 2. Relationship between carbon densities of typical vegetation and stand age. 

doi:1 0.1 371 /journal.pone.0082208.g002 



vegetation survey, organic carbon analysis data, and NFRI, the 
carbon densities of typical vegetations were fitted against age class: 



Using the carbon contents obtained through laboratory analysis, 
the SOC per hectare was calculated as [35] : 



--ax h 



(2) 



where x is the stand age, a and b are the parameters to be fitted. The 
fitted growth curves (Fig. 2) were used to estimate the carbon density 
of typical vegetations in the target year. The curves represent the 
broadleaved, coniferous, and economic forests, and the shrubs. 

The biomass densities of shrubs, herbs, and litters were obtained 
by multiplying biomass expansion factor (BEF) by dry plant weight 
per unit quadrat, and BEF was calculated as follows: 



BEF- 



10000m 2 



( (/»(/ sizc(m 2 ) 



(3) 



Soil bulk density was calculated as follows [35-36] : 

Bulk density{g-cm ) = 

Oven dry mass(g cm~ ) 

, / Mass of coarse fragments(g) 
Core volumeicm 5 ) — — : — 



\Density of rock fragments(g-cm 3 ) 



_(4) 



SOC(Mg-ha~ 1 ) = 

[(«>;7 bulk density(g-cm~ 3 )xsoil depth(cm)xcarbon content(g-kg 3 ))] 
x 100 (5) 



The average SOC changing rate from afforestation was 
calculated as a time-weighted mean value as follows [10], [37]: 



C time — weighted ^ ^ C se q-J ^ ^ age t 



(6) 



where C seqi and age, are the changes in SOC density and year 
respectively after afforestation; n is number of observations. 

An area-weighted total carbon density was then calculated for 
the broadleaved, coniferous, and economic forests and the shrubs 
separately. Carbon pools were estimated as the sum of total carbon 
density multiplied by forest category and stand size-class, and the 
corresponding area [38]. The data for the four forest types were 
combined to estimate total carbon storages of plantation 
ecosystems in sand source areas of north Beijing. 
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Data analyses 

One-way analysis of variance (ANOVA) was used to detect 
statistical differences between carbon stocks or carbon densities 
among forest types or regions. Multiple comparisons were made 
using Duncan's multi-range test in cases the data were in a 
significandy non-normal distribution. Unless otherwise indicated, 
an alpha of 0.05 was considered significant. Analyses were 
performed using SPSS 13.0 on Windows. 

Results 

Area changes 

The afforestation area in sand source areas of north Beijing was 
2.3 million ha (Table 2), or 1.16 percent of China's forest area (195 
million ha, 2004-2008 NFRI). The plantation areas in 
WRPZYMHR, DLZAPR, and OSLZ were 1.1, 0.4, and 0.8 
million ha, respectively (Table 2). The areas of broadleaved, 
coniferous, and economic forests, and shrubs were 0.5, 0.8, 0.4, 
and 0.6 million ha, respectively (Table 2). 

Carbon storage in plantation ecosystems in sand source 
areas 

Plantation ecosystems in sand source areas of north Beijing 
contained about 173.02 Tg C (Table 2), or 2.22 percent of the 
total carbon stored in China's forests (7.81 Pg C, [39]). About 
9.99 Tg C was stored in vegetations, and the rest was stored in 
soils (163.03 Tg C) (Table 2). 

The average carbon storage was 93.12 Mg ha -1 in plantation 
ecosystems (Tables 3-4), and varied among the components. Trees 



including roots contained 8.42 Mg ha of carbon, or 9.04 
percent of total forest ecosystem carbon (Table 3). The shrub-herb 
layer and the litter layer contained 5.55 and 2.05 Mg ha -1 of 
carbon respectively, or 5.96 and 2.20 percent of total forest 
ecosystem carbon, respectively (Table 3). The largest proportion of 
carbon was in the soil, 77.10 Mg ha , or 82.80 percent of the 
total forest ecosystem carbon (Table 4). Soil carbon storage of the 
four forest stands decreased with soil depth, but it was mainly 
distributed in 0-40 cm depth of surface soil, accounting for 50.46 
percent of the total forest ecosystem carbon (Table 4). From the 
aspect of spatial distribution, the average forest ecosystem carbon 
storage decreased in the order: soil > tree > shrub-herb > litter. 
Soil was the main carbon pool for the four forest types. 

Carbon storages among regions 

Forest vegetation carbon storage was not distributed evenly 
across regions. Of the total 173.02 Tg C in plantation ecosystems, 
89.33 Tg C or 5 1 .63 percent was distributed in WRPZYMHR, far 
more than in DLZAPR (27.32 Tg C or 15.79 percent) and OSLZ 
(56.37 Tg C or 32.58 percent) (Table 2, Fig. 3A). 

The carbon quantity varied considerably among regions, with 
94.19 Mg ha" 1 in WRPZYMHR, 65.72 Mg ha" 1 in DLZAPR, 
and 102.38 Mg ha" 1 in OSLZ (Tables 3-4, Fig. 3B). The average 
vegetation carbon storage in WRPZYMHR was 13.91 Mg ha , 
significandy higher than in DLZAPR ( 1 1 . 1 4 Mg ha" ') and OSLZ 
(7.29 Mg ha l ) (Table 3, Fig. 3B). The majority of carbon was 
stored in the soil pool, with 80.28, 54.58, and 95.09 Mg ha" 1 
stored in WRPZYMHR, DLZAPR, and OSLZ, respectively 



Table 2. Carbon storages of plantation ecosystems in 


sand source 


areas of North Beijing, China (mean ± 


S.E.). 








Area (million 
ha) 


Vegetation carbon Soil carbon storage 
storage (Tg C) (Tg C) 


Total carbon storage 
(Tg C) 


Water resources protection zone in Yanshan 
mountainous and hilly region 


Broadleaved 
forests 


0.1 


1.26±0.05 


10.66±1.01 


11.92±2.34 




Coniferous 
forests 


0.6 


3.08±0.10 


47.76±4.34 


50.84±6.54 




Economic 
forests 


0.2 


0.41 ±0.00 


12.88±1.11 


13.29±2.34 




Shrubs 


0.2 


0.67±0.00 


12.61 ±2.21 


13.28±2.21 




Total 


1.1 


5.41 ±0.45 


83.92±9.89 


89.33±9.89 


Desertificated land zone in agro-pasture region 


Broadleaved 
forests 


0.0 


0.51 ±0.00 


4.22±0.98 


4.73 ±1.23 




Coniferous 
forests 


0.1 


0.61 ±0.00 


10.84±1.23 


11.44±3.32 




Economic 
forests 


0.0 


0.03 ±0.00 


1.56 ±0.00 


1.59 ±0.00 




Shrubs 


0.2 


0.51 ±0.00 


9.05 ±0.03 


9.56 ±0.67 




Total 


0.4 


1.65±0.03 


25.67±5.65 


27.32±7.88 


Otingdag sandy land zone 


Broadleaved 
forests 


0.3 


1.75±0.02 


23.91 ±3.45 


25.66±4.54 




Coniferous 
forests 


0.1 


0.41 ±0.00 


9.00 ±0.02 


9.41 ±2.32 




Economic 
forests 


0.2 


0.54±0.00 


12.55±0.34 


13.09±3.21 




Shrubs 


0.2 


0.24±0.00 


7.97 ±0.45 


8.21 ±1.21 




Total 


0.8 


2.93±0.12 


53.44±7.67 


56.37±9.98 


Total 




2.3 


9.99±0.34 


163.03±14.23 


173.02± 13.21 



doi:1 0.1 371 /journal.pone.0082208.t002 
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(Table 4, Fig. 3B). The highest average carbon in soils was in 
OSLZ. 

Carbon storages among forest types 

Total carbon storage and the relative storage by forest 
ecosystem components were both significantly different among 
the four forest types (Table 3). Biomass density ranged from 52.98 
to 278.19 (are a- weighted mean 108.93) Mg ha 1 in broadleaved 
forests; from 49.72 to 229.81 (109.34) Mg ha" 1 in coniferous 
forests; from 35.03 to 1 10.43 (66.66) Mg ha 1 in economic forests; 
and from 22.38 to 109.61 (64.79) Mgha" 1 in shrubs (Fig. 4A). For 
the total 173.02 Tg C in plantation ecosystems, broadleaved, 
coniferous, and economic forests and shrubs sequestrated 42.30, 
71.69, 27.98, and 31.05 Tg C, respectively (Table 2, Fig. 4B). Both 
L. gmelini and H. rhamnoides forests contained high amounts of 
carbon, with 23.44 and 21.98 Mg ha 1 carbon stored in 
vegetation respectively, and 165.49 and 78.74 Mg ha 1 in 
vegetation/soil respectively (Table 3-4). Both U. pumila and A. 
sibirica forests contained low amounts of carbon, with 4.21 and 
6.23 Mg ha 1 stored in vegetation respectively, and 83.86 and 
71.46 Mg ha 1 respectively in aboveground/belowground parts 
(Table 3-4). 

Discussion 

From 2001 to 2010, the carbon storage in the 2.3 million ha of 
plantation ecosystems in sand source areas of north Beijing was 
173.02 Tg, or 2.22 percent of all the carbon stored in China's 
forests. Plantation ecosystems in sand source areas play an 
important role in the carbon cycle and carbon balance of 
terrestrial ecosystems in China. The increased forest area and 
the plantation growth suggested that successful afforestation 
implementation and consecutive management programs over the 
past decade allow plantations to contribute greatly to carbon 
sequestration [40]. Results from China and other countries 
indicated that plantations made major contributions to the 
formation of state-level carbon sink [39-41]. China's plantation 
area is 62 million ha according to 2004—2008 NFRI, and forest 
ecosystems have turned into carbon sinks because of the state-level 
afforestation and reforestation activities during the past few years 
[39]. Plantations in China and Japan made major contributions to 
state-level carbon sink [41-42]. The total forest biomass carbon 
stocks in Korea increasing from 20.57 Tg in 1954 to 239.85 Tg in 
2007 were mainly attributed to the successfully implemented 
afforestation and subsequent management practices [40]. Mean- 
while, some phenomena were also found in the course of project 
implementation and management in this research. For example, in 
our study areas, economic forests developed higher levels of 
vegetation cover than other afforestation species over the past 
decade. The undergrowth herbaceous vegetations were direcdy 
removed to accelerate tree growth, which is the initiative for forest 
restoration [32]. Species diversity which plays a vital role in 
enhancing carbon storage and overall productivity of plantations 
[43] was seriously diminished because of management activities 
and artificial tending. Extensive afforestation efforts in arid and 
semiarid regions have caused environmental deterioration, with 
deteriorating ecosystems and more frequent wind erosion [44]. If 
species selection or early land-shaping and artificial tending is 
inappropriate, even traditional forms of afforestation can be risky, 
and plantations may fail [45—46]. Therefore, preserving or 
restoring the existing ecosystems by conserving and actively 
managing the large areas of existing secondary forests may be a 
primary goal of carbon sequestration efforts [46-47], and many 
degraded ecosystems demonstrate a remarkable ability to recover 
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Figure 3. Carbon storages (A) and carbon densities (B) in plantation ecosystems among regions. I: Water resources protection zone in 
Yanshan mountainous and hilly region; II: Desertificated land zone in agro-pasture region; III: Otingdag sandy land zone. 
doi:1 0.1 371 /journal.pone.0082208.g003 



through natural processes rather than afforestation [48] . However, 
the natural recovering ability of severely degraded forests cannot 
be assured, and recovery becomes more difficult when an 
ecosystem has crossed an ecological threshold to a new steady 
state [49]. Moreover, various afforestation methodologies are 
available for currendy deforested areas [46]. Future afforestation 
designs should target more at selection of suitable tree species, 
management of age structures, and implementation of good 
silvicultural practices for improving biodiversity, ecological func- 
tioning, and human livelihoods [50-51]. 

Forest carbon storage in a specific region is affected by many 
factors such as climate, solar radiation, forest area, forest ages, 
species composition, land-use history, disturbance, and field and 
soil characteristics [51-52]. This study shows that coniferous 
forests have statistically higher total carbon density than other 
forest types. Among all coniferous forests, L. gmelini forest contains 
the highest average carbon because of the large carbon quantity 
stored in trees, which is consistent with other studies [53-54], so it 
may be concluded that coniferous forests are the plant species 



more suitable for the local environment. Economic forests such as 
A. sibirica contain a lower understory carbon because the 
undergrowth herbaceous vegetation (i.e. grasses, forbs, and herbs) 
was manually removed to promote tree growth by alleviating 
competition for nutrient and moisture. Among broadleaved 
forests, the average carbon density of U. pumila is lower because 
the slow-growing younger plantations support lower vegetation 
biomass. Among shrubs, H. rhamnoides contains higher carbon 
because its rapid root propagation supports higher vegetation 
densities. The differences in forest area, forest species, and forest 
structures led to different ecosystem carbon sequestrations [32], 
[51]. Increasing forest area and vegetation carbon density through 
application of different tree species and age structures can be a 
practical and effective strategy to improve vegetation carbon 
sequestration [51], [55]. 

The carbon quantity varied considerably among regions, as the 
average vegetation carbon storages in WRPZYMHR, DLZAPR, 
and OSLZ were 13.91, 11.14, and 7.29 Mg ha -1 respectively, and 
the average soil carbon storages were 80.28, 54.58, and 95.09 Mg 




Figure 4. Carbon densities (A) and carbon storages (B) in plantation ecosystems among forest types. BF: Broadleaved forests; CF: 
Coniferous forests; EF: Economic forests; S: Shrubs. 
doi:1 0.1 371 /journal.pone.0082208.g004 
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ha , respectively (Tables 3-4, Fig 3B), Zhou et al. [56] and 
Huang [57] also found that carbon density in the forest ecosystem 
was spatially differentiated with the zonalities of longitude, latitude 
and altitude, probably because of the differences in regional 
climates and soil conditions. The results also suggested an 
important direction for design of reforestation programs and for 
enhancement of vegetation carbon sequestration potentials in sand 
source areas of north Beijing [51]: target more at the selection of 
suitable tree species according to local climate and soil conditions. 

The average carbon density in plantation ecosystems is 
93.12 Mg ha -1 , and carbon storage in forest vegetation is 
25.06 Mg ha" , which are generally lower compared with other 
studies in China [58-61]. For instance, the average carbon density 
in forest vegetation in China was determined to be 42-50 Mg 
ha" 1 (15-32 Mg ha" 1 in planted forests) [58], 35-39 Mg 
ha _1 [61], or 6.47-1 18.14 Mg ha" 1 [60]. The main cause for 
the differences might be that the plantations in sand source areas 
of north Beijing are younger. It also suggests that the forests there 
have large potential for carbon sequestration. In this study, forest 
carbon storage includes the carbon stored in trees, understory, 
forest floor, and top soils (i.e. 0-100 cm). Trees including roots 
contained 8.42 Mg ha 1 of carbon, while understory, litter layer, 
and the 0-100 cm topsoil contained 5.55, 2.05, and 77.10 Mg 
ha" 1 , respectively, which are generally lower than that of other 
forested areas in China. For instance, the average carbon densities 
of vegetation, litter, and soil in Chinese forest ecosystems were 
estimated to be 57.07, 8.21, and 193.55 Mg ha" 1 , respectively 
[56]. Lower carbon accumulation in the studied forest and soil 
might result from the singleness of afforestation species, younger 
forest age, lower biomass, poorer forest structure and function, and 
shorter restoration time [32], [62]. Land management practices, 
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including species selection, land-shaping patterns, and artificial 
tending patterns should be adjusted in design of afforestation 
programs to increase vegetation carbon storage. Soil carbon 
storage distributed in the top 0-40 cm layer accounts for 50.46 
percent of the total carbon storage within the 0-100 cm of soil, 
which is consistent with other studies [32], [57], [63]. Reportedly, 
more than 60 percent [32], or 35-80 percent of carbon in the 
depth of 0-100 cm was stored in the top 0 — 40 cm in tropical and 
subtropical soils [64] . Therefore, protection of organic carbon in 
the surface soils of planted forests is important for improving 
carbon sequestration [32]. 

Conclusions 

Younger plantations have large potential for enhancing 
terrestrial ecosystem carbon sequestration. In future afforestation, 
it will be more effective to focus on increasing forest area and 
vegetation carbon density through selection of appropriate tree 
species and stand structure according to local climate and soil 
conditions, and through application of proper forest management 
including land-shaping, artificial tending and fencing plantations. 
It will be also important to protect the organic carbon in surface 
soils during forest management. 
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